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Besides perturbation theory, which requires the knowledge of the exact unperturbed solution, 
variational techniques represent the main tool for any investigation of the eigenvalue problem 
of some semibounded operator H in quantum theory. For a reasonable choice of the employed 
trial subspace of the domain of H, the lowest eigenvalues of H can be located with acceptable 
precision whereas the trial-subspace vectors corresponding to these eigenvalues approximate, 
in general, the exact eigenstates of H with much less accuracy. Accordingly, various measures 
for the accuracy of approximate eigenstates derived by variational techniques are scrutinized. 
In particular, the matrix elements of the commutator of the operator H and (suitably chosen) 
different operators with respect to degenerate approximate eigenstates of H obtained by the 
variational methods are proposed as new criteria for the accuracy of variational eigenstates. 
These considerations are applied to that Hamiltonian the eigenvalue problem of which defines 
the spinless Salpeter equation. This bound-state wave equation may be regarded as the most 
straightforward relativistic generalization of the usual nonrelativistic Schrodinger formalism, 
and is frequently used to describe, e.g., spin-averaged mass spectra of bound states of quarks. 
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1 Motivation 

A central element in quantum theory is the solution of eigenvalue problems. However, 
usually there is no suitable exact solution to which perturbation theory can be applied. 

A very efficient way to locate the discrete spectrum of some self-adjoint operator H 
bounded from below is provided by the famous Rayleigh-Ritz variational technique |l|]: 
If the eigenvalues Ek, k = 0,1, of H are ordered according to Eq < Ei < E2 < . . . , 
the first d of them are bounded from above by the d eigenvalues Ek, k = 0,1, . . . ,d — l, 
(ordered by Eq < Ei < . . . < E^-i) of that operator which is obtained by restricting H 
to some d-dimensional subspace of the domain of H, i.e., Ek < Ek, k = 0,1, . . . ,d — l. 
Unless the eigenstate \xk) corresponding to some eigenvalue Ek is already an element of 
this trial space — in which case the upper bound Ek becomes identical to the eigenvalue 
Ek — , enlarging the dimension d of the trial space will, in general, improve the obtained 
upper bound; at least, it can't make things worse. However, it is not straightforward to 
quantify how close approximate eigenstates resulting from the variational method and 
exact eigenstates are.QThus, we embark upon a systematic study of the accuracy of the 
variationally determined eigenstates of H and suitable measures to judge their quahty. 



2 Measures of Quality for Variational Trial States 

Consider some self-adjoint operator H, W = H , assumed to be bounded from below. 
Suppressing, for the moment, the index A; = 0, 1, 2, . . . discriminating between different 
solutions, let the eigenvalue equation for H, 

H\x)=E\x), (1) 

be solved by some (generic) eigenvector \x) corresponding to some (real) eigenvalue E, 
to be extracted according to 

j,^ {x\H\x) 

ixlx) 

The Rayleigh-Ritz variational technique yields an upper bound E on this eigenvalue E 
as well as, by diagonalization of the relevant characteristic equation, the corresponding 
vector \(p) in the rf-dimensional trial space. There exist several (potentially meaningful) 
measures of the quality of this trial state \ip) which immediately come to one's mind: 

1. The trial state \ip) is supposed to represent — to a certain degree of accuracy — the 
approximate solution of the eigenvalue problem defined in Eq. ([^) . Consequently, 
a ffist indicator for the resemblance of \(f) with the exact eigenstate \x) would be 
the distance between the expectation value of the operator H with respect to the 
trial state \(p), i.e., between the obtained upper bound 

and the exact eigenvalue E. However, the precise location of the exact eigenvalue 
E is usually not known. Lower bounds on eigenvalues of self-adjoint operators are 
much harder to find than upper bounds. The practical use of Temple's inequality 
[0 is, unfortunately, limited. The "local-energy" theorem applies, in general, 
only to the ground state. 

^ Recent investigations of the production and decay of heavy quarkonium systems have stimulated 
renewed interest in the exact value of the wave function at the origin of the two-quark bound state 12] . 
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2. The natural measure for the resemblance of the Hilbert-space vectors \ip) and \x) 
under consideration is the overlap 



of the trial state Iv?) with the eigenstate If the generic eigenstate |x) refers to 
the ground state |xo) of the operator H, the deviation of this overlap from unity 
is, according to the (rather efficient) Eckart criterion [Q, bounded from above by 

^ -T-T-\ ^0 = • (3) 



A comprehensive discussion of both lower and upper bounds on the overlap 
S as criteria for the accuracy of approximate wave functions may be found 
in Ref. f^. As is evident even from the simple Eckart criterion all these 
bounds require, in general, some additional knowledge, such as the location of 
the eigenvalues i^^, A; = 1, 2, . . . , of the operator H or matrix elements of higher 
powers if", n > 2, of H, for instance, certain moments 

(y^l^-ly^) ^>2 

of H, or, at least, appropriate bounds on these quantities 0. Rather frequently, 
however, the required information is either not available at all or hard to obtain. 

Consider the commutator [G, H] of the operator H under consideration with any 
other operator G, where the domain of G is assumed to contain the domain of H. 
Then the expectation value of this commutator with respect to a given eigenstate 
\x) oi H (or, more generally, the matrix element of this commutator with respect 
to an arbitrary pair of degenerate eigenstates \xi) and \xj) of H, i.e., eigenstates 
which satisfy Ei = Ej) trivially vanishes: 

{x\[G,H]\x) = 0. (4) 

(Under the name "hypervirial theorems" 0, this relation has been generalized, 
for Hamiltonians H which involve nonrelativistic kinematics, to nonnormalizable 
vectors \x), that is, to vectors which are not elements of some Hilbert space [0.) 
Hence, choosing different operators G generates a whole class of operators [G, H] 
each of which may serve to test the quality of a given trial state \ip) by evaluating 
how close the expectation value (v?|[G, H]\ip) of this operator with respect to \ip) 
comes to zero: 

{cp\[G,H]\^)^0. 

This expectation value vanishes, of course, also if, by accident, the state \ip) is an 
eigenstate of G. However, for a given operator G, after having determined \ip), it 
is straightforward to check for this circumstance, for instance, by inspecting the 
variance Aq of the operator G with respect to the state 1^9), defined, as usual, by 

f {v\G\^) y ^ 
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which clearly vanishes if the state \(p) under consideration is an eigenstate of G. 
Moreover, it goes without saying that an expectation value (v^l [G, H]\ip) vanishes 
also if the state \ip) is an eigenstate of the commutator [G, H] with eigenvalue 0, 
or even if the state defined by [G, H]\({)) proves to be orthogonal to the state \({)). 

For any self-adjoint operator G, i.e., G^ = G, this commutator is anti-Hermitean, 
which obviously suggests to define a self-adjoint operator G = G"'' (on the domain 
of the operator H) by 

[G,H] =: iC . 

If, for example, G is chosen to be the (symmetrized and self-adjoint) generator of 
dilations, 

G = ^ (x ■ p + p ■ x) , (5) 

the relation (^) is precisely the "master virial theorem" introduced in Ref. for 
a systematic study of (relativistic) virial theorems 0. In this case, for operators 
H of the form of some typical Hamiltonian consisting of a momentum-dependent 
kinetic-energy operator, T(p), and a coordinate-dependent interaction-potential 
operator, V{x.), that is, 

H = T(p) + y(x) , 
the operator G becomes the "virial operator" 

The point spectrum (i.e., the set of all eigenvalues) of the dilation generator (^ is 
empty; in other words, the dilation generator has no Hilbert-space eigenvectors. 



In any case, however, one should always keep in mind one (trivial) fact [10|: Depending 
on one's particular choice of the employed trial subspace, it may happen that one trial 
state l^'*-^'') yields a better approximation to the exact eigenstate \x), as judged by their 
overlap 

^(^a)|^a))(x|x) ' 

but a worse upper bound E^^^ on the corresponding eigenvalue E while some other trial 
state Iv?*-^-*) yields the better upper bound E'^'^^ < E^^^ on the eigenvalue E but a worse 
approximation to the exact eigenstate which fact would be betrayed by the overlap 



3 Prototype of Relativistic Wave Equations: The 
Spinless Salpeter Equation 

Let us apply the above general considerations to the prototype of all (semi-) relativistic 
bound-state equations, the "spinless Salpeter equation," defined by a — by assumption, 
self-adjoint — Hamiltonian H (in one-particle form, which encompasses the equal-mass 
two-particle case, too |jll|, |12|, |13|, [14|): 



H = T + V , (7) 
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where T is the "square-root" operator of the relativistic kinetic energy of some particle 
of mass m and momentum p, 



and V = V^(x) is an arbitrary, coordinate-dependent, static interaction potential. The 
spinless Salpeter equation is then just the eigenvalue equation for this Hamiltonian H, 

H\xk) = Ek\xk) , A; = 0,1,2,... , 

for the complete set of Hilbert-space eigenvectors \xk) of corresponding to its energy 
eigenvalues 

^ ^ {Xk\H\xk) 
{Xk\Xk) 

Analytic upper bounds on these eigenvalues have been given p, ii2|, m, m m m 



For the Coulomb potential, the local-energy theorem has been successfully applied [p!8[] . 

For the sake of comparison, we focus our interest to central potentials V{yi) = V{r), 
r = |x|. Furthermore, in order to facilitate the numerical treatment of the problem, we 
consider the harmonic-oscillator potential 

V{r) = ar'^, a>0. (9) 

The reason for this particular choice is the following: In momentum space, the operator 
is represented by the Laplacian with respect to the momentum p, — > — Ap, while 
the kinetic energy T, nonlocal in configuration space, is represented by a multiplication 
operator. Consequently, exactly for a harmonic-oscillator potential the semirelativistic 
Hamiltonian H in its momentum-space representation is equivalent to a nonrelativistic 
Hamiltonian with some (effective) interaction potential reminiscent of the square root: 



= -aAp + ^p2 + m2 . (10) 

The solutions of the corresponding eigenvalue equation may then be found with one of 
the numerous procedures designed for the treatment of the nonrelativistic Schrodinger 
equation. 

For the harmonic-oscillator potential, it is comparatively easy to get a first idea of 
the approximate location of the energy levels Ej. by entirely analytical considerations: 

• On the one hand, since 

p2 

r(p) < m + — , 
2m 

every energy eigenvalue of H is bounded from above by the energy eigenvalue 
Ek,T<}R of the nonrelativistic counterpart 

i^NR = m + h V 

2m 

of H: 

Ek < -Efc,NR ■ 

For the harmonic-oscillator potential (pi), these nonrelativistic energy levels read 




E^n = m+\l'^ (2n, + i + ^y (11) 

expressed in terms of radial quantum number n,. = 0,1,2, .. . and orbital angular 
momentum i = 0,1,2, ... . An (easily) improved bound may be found in App. ^ 
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• On the other hand, since 

Tip) > IpI , 

every energy eigenvalue Ek of H is bounded from below by the energy eigenvalue 
Ek{m = 0) of the Hamiltonian 

H{m = 0) = \p\ + V 

corresponding to vanishing particle mass m: 

Ek > Ek{m = 0) . 

For the harmonic-oscillator potential (j^), by remembering the momentum-space 
representation ([10|) of H all bound states of vanishing orbital angular momentum 
i ( "S waves" ) may be obtained as solutions of a differential equation of the form 

the resulting energy levels read 

E{m = 0) = -a^^^zo , 

where Zq denote the zeros of the Airy function Ai{z), which solves this differential 
equation: zo = -2.33810 . . . , -4.08794 . . . , -5.52055 ....... Energy levels 

of bound states of nonvanishing orbital angular momentum i ("P, D, . . . waves") 
must be determined numerically anyway. Also this bound is improved in App. ^ 

For homogeneous functions, the radial derivatives entering in the virial operator (|^) 
simply probe the respective degree of homogeneity of these functions. In particular, for 
the harmonic-oscillator potential @ one finds 

X - — V(x) =r—V{r) = 2ar^ = 2V{r) . 
ax or 

With the expression for C, introducing the (only momentum- dependent) operator 

F.T(p) + lp.|-T(p). 
the Hamiltonian with harmonic-oscillator potential thus may be cast into the form 

H = r(p) + v(x) 

. T(p) + ix.Anx) 

- r(p) + 5P-|;^c)-f 

2 

This observation allows to express any matrix element of the operator C (in particular, 
its expectation value {lp\C\lp) with respect to the trial state \ip) under investigation) as 
the corresponding matrix element of the difference of Hamiltonian H and operator F: 

{cp\C\ip) =2{ip\F-H\^)^2 {{ip\F\ip) - E {iflif)) . 

For the relativistic expression (P) of the kinetic energy T, the above operator F reads 

^ 3 p2 + 2 



2 Jp^ + m2 
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4 Definition of the "Laguerre" Trial Space 

As far as the achieved accuracy of the solutions obtained is concerned, the most crucial 
step in all variational games of the Rayleigh-Ritz kind is, for a given operator H under 
consideration, a reasonable definition of the adopted trial subspace of the domain of H . 

For spherically symmetric potentials V{r), a very popular choice for the basis states 
which span the trial space required for the application of the variational technique are 
"Laguerre" trial states, given in configuration-space representation by |^, p!4| 



where L^^'' {x) are the generalized Laguerre polynomials (for parameter 7) |]T9[ , defined 
by the power series 

4-w-|:M)'(r-])^ 

and normalized, with the weight function x"^ exp(— a;), according to 

r(7 + A; + 1 



°° i ' -|- fc -|- i 
j dxx^exp(-x) L^k\x) lJ.T^(x) = ^ 5kk' 



and y^rn{^) are the spherical harmonics for angular momentum ^ and its projection m, 
depending on the solid angle VL and orthonormalized according to 

The trial functions (|12|) involve two variational parameters, /x (with dimension of mass) 
and (3 (dimensionless), which, by the requirement of normalizability of these functions, 
are subject to the constraints /i > and 2/3 > —1. 

One of the advantages of the trial function (|12]) is the easy availability of an analytic 
expression for the corresponding momentum-space representation of these trial states, 
obtained by Fourier transformation: 



^fe/m(p) 



(2 /i)2^+2/5+i k\ (-i)^ IpI^ ^ (-1)* ( k + 2^ + 2^3 



i-iY 


|P 




21+1/2 Y ( 





i=0 



t! I k-t 



T{2i + (3 + t + 2){2fiy ^f2e + f3 + t + 2 + 3 

-r , — - — + 



(p2 + /i2){2£+/3+t+2)/2 \^ 2 ' 2 ' 2'p2+/i2; 

with the hypergeometric series F, defined by 

T{w) ^ T{u + n)T{v + n) z"" 

F{u, v; w; z) = , , , , }^ ■ — ^ • 

r(M)r(t;)^Q T{w + n) nl 

For the present investigation, we too employ the "Laguerre" trial states defined by 
Eq. (|T2p, with, for both definiteness and ease of calculation, the variational parameters 
/i and P kept fixed to the values n = m and /? = 

^ Within some specific nonrelativistic quark potential model, the significance of the wave function 
at the origin obtained variationally from superpositions of Gaussian trial functions has been discussed 
numerically [El|. 
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Of course, the trivial fact that the (Hilbert-space) states defined by the "Laguerre" 
trial functions ( |T2l) cannot be eigenvectors of the dilation generator (|]) is reflected by 
the nonvanishing numerical value of the variance Ac of the operator G in these states: 
For instance, for the lowest-dimensional possibility, d = 1, realized by the trial function 

^o,oo(x) = ]j^exp{-nr) , 

this variance is exactly = 3/4. For trial-space dimension d = 25, this quantity takes 
in the state corresponding to the lowest eigenvalue E the numerical value A^ = 1.3779. 

Up to a dimension d = 4 of the trial space, the matrix elements of the Hamiltonian 
H with respect to the chosen trial states are, at least in principle, accessible by entirely 
algebraic manipulations. In the simplest conceivable case, realized by the choice d = 1, 
one finds with the help of the explicit analytic expressions for the matrix elements of H 
derived in Ref. |jT6|, for instance, for the ground-state energy Eq the naive upper bound 

- 64 m 3 a 
-^0 = ^ 2 ■ 

5 Rates of Convergence of the Quality Measures 

Now, let us observe our variational eigenstates, \(p), approaching the exact eigenstates, 
\x), for increasing dimension d of the employed trial space, by comparing the behaviour 
of the various measures for the accuracy of approximate eigenstates discussed in Sec. ^ 

Without doubt, the only genuine "point of reference" of any variational solution to 
an eigenvalue problem is the corresponding exact solution. The exact solution sought is 
computed here with the help of the numerical integration procedure developed for the 
solution of the nonrelativistic Schrodinger equation in Ref. (Since the problematic 
term in our Hamiltonian is the kinetic-energy operator the solution is constructed 
in momentum space,[| which, in addition, allows to take advantage of the operator F.) 

Table ^confronts, for the ground state and the lowest radial and orbital excitations, 
the approximate solutions as calculated with the help of the Rayleigh-Ritz variational 
technique for "Laguerre" trial subspaces of the domain of H of increasing dimension d 
with the exact solutions of the eigenvalue problem for the semirelativistic Hamiltonian 

with a central interaction potential of the harmonic-oscillator form (|^). First of all, 
as discussed in Sec. the exact position of any eigenvalue E of our Hamiltonian H is 
confined to a range defined by the nonrelativistic upper bound i?NR and the zero-mass 
lower bound E{m = 0) on this energy eigenvalue E. There are several quantities which 
may participate in a competition for "the best or most reasonable measure of quality:" 

1. The relative error e of every (Rayleigh-Ritz) upper bound E on the exact energy 
eigenvalue E, 

E-E 

e = — ^ , (13) 
is, by definition, always nonnegative, i.e., 

e>0 . 



The outcome of any numerical integration procedure is a solution given only at discrete points in 
the respective representation space. The reconstruction of the entire function, with the help of some 
interpolation procedure, in form of an interpolating function introduces, clearly, an additional error. 
We have, of course, checked that this error is negligible. 
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Table 1 : Characterization of the quahty of the variational solution of the eigenvalue problem 
of the semirelativistic Hamiltonian H = ■\/p^ + rr?+y(r) with harmonic-oscillator potential 
V{r) = ar'^, for states of radial quantum number = 0, 1, 2 and orbital angular momentum 
£ = 0,1,2 (called IS, 2S, 3S, IP, and ID in usual spectroscopic notation), obtained with the 
help of "Laguerre" trial states spanning trial spaces of increasing dimension d = 1,2, 3, 10, 25, 
by: the nonrelativistic upper bound -Enr and zero-mass lower bound E{m = 0) on the energy, 
the (numerically computed) "exact" energy E, the variational upper bound E on this energy, 
the relative error e of the upper bound, the deviation from unity, a, of the overlap squared of 
exact and variational eigenstates, the (appropriately normalized) expectation values of the 
virial operator C, and the (normalized) maximum local difference lv of the momentum-space 
representations of exact and variational eigenstates. The physical parameters are fixed to the 
values m = 2 GeV for the particle mass and a = 2 GeV'^ for the harmonic-oscillator coupling. 
A simple entry "0" indicates that the numerical value is closer to than the rounding error. 
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e [of. Eq. (y)] 


1 


0.1023 






0.3280 


0.6583 




2 


0.0268 


0.4002 




0.0694 


0.2180 




3 


0.0266 


0.1059 


0.9294 


0.0178 


0.0722 




10 


0.0001 


0.0050 


0.0212 


0.0005 


0.0008 




25 








0.0001 








a [of. Eq. (Ill)] 


1 


0.09618 






0.36144 


0.65587 




2 


0.02375 


0.43693 




0.09001 


0.34398 




3 


0.02280 


0.13878 


0.83034 


0.01918 


0.12705 




10 


0.00003 


0.00560 


0.03727 


0.00040 


0.00061 




25 








0.00008 








y [of. Eq. (©] 


1 


-0.6120 






-0.8328 


-0.9074 




2 


+0.0308 


-0.8666 




-0.5103 


-0.7483 




3 


-0.0187 


-0.4577 


-0.9525 


-0.1275 


-0.5197 




10 


-0.0034 


-0.0167 


-0.1895 


-0.0016 


-0.0175 




25 





-0.0001 


+0.0001 








u [cf. Eq. (0)] 


1 


+0.9277 






+0.7541 


+1.0578 




2 


-0.00754 


+2.4577 




+0.3598 


+0.7262 




3 


-0.01049 


+0.2265 


+4.3252 


+0.1146 


+0.4209 




10 


-0.00867 


+0.0587 


+0.2592 


+0.0183 


+0.0186 




25 


+0.00003 


-0.0017 


+0.0002 


+0.0004 


+0.0003 
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2. The deviation from unity, cr, of the modulus squared of the overlap S of exact and 
variational eigenstates defined in Eq. (H), 

a^l-\S\\ (14) 

is clearly confined to the range 

< (7 < 1 ; 

the lower bound applies when the variational eigenstate \(f) agrees with the exact 
eigenstate \x), while the upper bound applies when the variational eigenstate \ip) 
is orthogonal to the exact eigenstate Just for comparison, the Eckart bounds 
(H) on a, unfortunately available only for the ground state, are given in Table 
In our case, these bounds turn out to be always larger than the actual numerical 
values of the ground-state overlap by at least a factor 2. 

Table 2: Eckart bounds on the difference to unity of the overlap squared a of Eq. ( p^ ) (which 
apply only to the ground state, identified by vanishing radial and orbital-angular-momentum 
quantum numbers = i = 0) for increasing dimension d = 1,2,3, 10, 25 of the trial subspace 









-Eo 




-Eo 



1 0.19903 

2 0.05222 

3 0.05183 
10 0.00019 
25 0.00001 



3. The use of the expectation values of the commutators [G, H] with respect to the 
variational eigenstates \(p) is illustrated for the particular example of the dilation 
generator G defined in Eq. (|^), by considering (suitably normalized) expectation 
values {(f\G\(p) of the virial operator G given in Eq. (^: 

^ {^\G\cp) _ (¥^|p-iT(p)|v.) 

4. Finally, the normalized maximum difference of the normalized momentum-space 
representations <^(p) and x(p) of variational eigenstate \ip) and exact eigenstate 
\x), respectively, i.e., the maximum pointwise relative error in momentum space, 

max[(^(p) - x(p)] 

^ = r-- (16) 

maxx(p) 

is listed. 

Note that the only measure for the accuracy of approximate eigenstates \ip) which does 
not require any information other than the one provided by the variational technique 
is u, the (normalized) expectation values of the commutator [G, H] with respect to \ip). 
Inspection of Table |I] reveals that u represents indeed a sensitive measure of quality: for 
increasing trial-space dimension d it converges to zero at roughly the same rate as both 
energy and overlap error, e and a, but makes more sense than a pointwise error like u. 
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6 Summary and Conclusions 

This work has been devoted to a systematic investigation of the quahty of solutions of 
the eigenvalue problem for some (semibounded, self-adjoint) operator H by variational 
methods. Clearly, even if the overlap of a variationally derived approximate eigenstate 
and the corresponding exact eigenstate is close to unity, this does not necessarily imply 
that expectation values of given operators (among which this operator H which defines 
the particular eigenvalue problem under consideration is just one) with respect to these 
two eigenstates are equal to the same degree of accuracy. For instance, a given operator 
might probe regions in representation space where the difference between approximate 
and exact solutions is of less importance for the overlap. (Only for functions f{H) of if, 
the difference between approximate and exact expectation values of / is, in general, of 
second order in the deviation of the approximate eigenstate from the exact eigenstate.) 
In view of this, various measures for the accuracy of approximate eigenstates have been 
discussed. Apart from the possible exception of some more or less pathological special 
cases (cf. our corresponding remarks in Sec. the vanishing of the expectation values 
of the commutator of H and any other well-defined operator, taken with respect to the 
approximate eigenstates, provides a useful set of criteria for estimating the significance 
of the variational solution. This has been illustrated by considering the commutator of 
the Hamiltonian of the spinless Salpeter equation — which represents the first step from 
the nonrelativistic Schrodinger formalism towards incorporation of relativistic effects 
and which is used for the semirelativistic description of the spin-averaged mass spectra 
of bound states of "constituent" quarks within the framework of potential models (for 



analyses of these attempts see, e.g., Refs. [23, 24, 25|) — and the generator of dilations 
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A The Relativistic Harmonic-Oscillator Problem 

In spite of the fact that it is somewhat off the main stream of the present investigation, 
let us demonstrate how, with only slightly more effort, the bounds on the energy levels 
of the relativistic harmonic oscillator derived in Sec. ^ may be considerably improved:^ 

• By considering the square of the (obviously self-adjoint) operator T — /i, where /i 
is an arbitrary real parameter (with the dimension of mass), one obtains a set of 



operator inequalities for the kinetic energy T given in Eq. (g) p6|, |I2|, |T^, |17 



T<t±II^ forall;.>0. 
- 2/i ^ 

These inequalities immediately translate into upper bounds on the energy levels 
Eki = 0, 1, 2, . . . , of the relativistic harmonic-oscillator problem, which involve 
the corresponding nonrelativistic energy eigenvalues -Ea;,nr recalled in Eq. ([TTD : 

Ek < -Efc,uppGr(/i) = — ^ ^ y ~ (-Efc,NR -m), n> . 



We thank the referee, presumably Andre Martin, for his "suggestion" to include this discussion. 
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• On the other hand, by taking into account, for instance, the obvious positivity of 
the square of the (self-adjoint) operator ^"^ — (1 — ^'^) for an arbitrary real 
parameter it is straightforward to convince oneself of the validity of the set of 
operator inequalities 

T > IpI + , 0<^<1. 

Remembering the energy eigenvalues Ek{m = 0) of the Hamiltonian H{m = 0) 
corresponding to vanishing particle mass m introduced in Sec. |^, these operator 
inequalities may be reformulated in a straightforward way as lower bounds on the 
energy levels E^, k = 0,1,2, of the relativistic harmonic-oscillator problem: 

Ek > ^fc,iower(0 = ^m + {l~ ey/^ Ek{m = 0) , < ^ < 1 • 

Numerical optimization of these bounds with respect to the parameters fi and ^ then 
yields the improved upper and lower bounds on the energies of the lowest-lying states 
of the relativistic harmonic oscillator listed in Table ^ 

Table 3: Improved upper and lower bounds on the lowest energy levels of the semirelativistic 
Hamiltonian H = \r^ + m^+V{r) with harmonic-oscillator potential V{r) = ar"^ , for states 
with radial quantum number n-^ = 0,1,2 and orbital angular momentum £ = 0, 1, 2 (denoted 
by IS, 2S, 3S, IP, and ID in usual spectroscopic notation), for a particle mass of m = 2 GeV 
and a harmonic-oscillator coupling strength of a = 2 GeV'^. 



Optimized Bounds 






State 








IS 


2S 


3S 


IP 


ID 







1 


2 








I 











1 


2 


inini;uppor(Ai) [GeV] 


3.89851 


5.99081 


7.80385 


4.98990 


5.99081 


max£;iowor(0 [GeV] 


3.75534 


5.68007 


7.36329 


4.85577 


5.86508 
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